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We study both the classical and the quantum target space of {j>, q) minimal string theory, 
using the FZZT brane as a probe. By thinking of the target space as the moduli space 
of FZZT branes, parametrized by the boundary cosmological constant x, we see that 
classically it consists of a Riemann surface Mp^q which is a p-sheeted cover of the complex 
X plane. However, we show using the dual matrix model that the exact quantum FZZT 
observables exhibit Stokes' phenomenon and are entire functions of x. Along the way 
we clarify some points about the semiclassical limit of D-brane correlation functions. The 
upshot is that nonperturbative effects modify the target space drastically, changing it from 
M.p,q to the complex x plane. To illustrate these ideas, we study in detail the example of 
(p, q) = (2, 1), which is dual to the Gaussian matrix model. Here we learn that the other 
sheets of the classical Riemann surface describe instantons in the effective theory on the 
brane. Finally, we discuss possible applications to black holes and the topological string. 
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1. Introduction 



Minimal string theories, or (p, q) minimal CFTs coupled to Liouville theory, are im- 
portant examples of tractable, exactly solvable models of quantum gravity. These models 
are interesting laboratories for the study of string theory because, despite their simplicity, 
they contain many of the features of critical string theory, including D-branes, holography 
and open/closed duality. First solved using the dual matrix model description [jl|-^ (for 
reviews, see e.g. PJTO[]), recent progress in the study of Liouville theory |1TT| - |T7[] has led to a 
greatly improved understanding of minimal string theory from the worldsheet perspective 

mm- 

One limitation of minimal string theory, however, has so far been the lack of a well- 
developed target space interpretation. In this paper, we will take the first steps towards 
a solution of this problem. Naively, the target space of minimal string theory is just the 
worldsheet Liouville field (p. However, it is common in string theory that different target 
spaces can have the same physics. In our case an equivalent description involves the free 
scalar field 4> related to (p by the non-local Backlund transformation (similar to T-duality). 
An important question which we will address is the distinction between the classical target 
space and the nonperturbative, quantum target space. We will see that they are quite 
different. 

Our point of view (which was used among other places in [^) is that a better effective 
description of target space can emerge out of the moduli space of D-branes. The advantage 
of this point of view is that it can capture all the nonperturbative corrections to the target 
space. We should point out, however, that different branes can lead to different target 
spaces. For example in compactification on a circle DO-branes probe the circle, while the 
Dl-branes probe the dual circle. 

Minimal string theories have D-branes (the FZZT branes) labelled by a continuous 
real parameter 

X = hb (1-1) 

(the boundary cosmological constant). We wish to interpret x as a target space coordinate, 
and the picture we have in mind is as follows. The minisuperspace wavefunction of the 
FZZT brane suggests that it is a D-brane in (p space stretching from (p = — oo and dissolving 
at ~ — ^ logx (where b = is the Liouville coupling constant). Therefore, the tip of 
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the FZZT brane at ~ — ^ logo; acts as a point-hke probe of the LiouviUe direction. It 
has the virtue of being able to penetrate into the strong-coupling region (j) +00, where 
one might expect there to be significant modifications to the classical target space. So 
in this description, target space is parametrized by the coordinate x. Large positive x 
corresponds to the weak-coupling region of the LiouviUe direction, while x of order one 
corresponds to strong coupling. 

To see how the worldsheet dynamics of the LiouviUe field modifies the naive target 
space, it is useful to analytically continue x to complex values. In the semiclassical approx- 
imation the D-branes are not single valued as a function of x and are labelled by a point 
in a finite multiple cover of the x-plane. This multiple cover corresponds to a Riemann 
surface M. which can be described as follows. In terms of the disk amplitude $ we define 

y = 9.$ (1.2) 

or equivalent ly 



X 



$ = y y{x')dx (1.3) 
Then, x and y satisfy an algebraic equation 

F{x,y) = Q (1.4) 

which describes the Riemann surface M. in C^. The parameters of the polynomial F{x, y) 
which determine the complex structure of M. depend on the parameters of the minimal 
string. For large x and y the equation F{x, y) = becomes of the form x*^ ~ y^ for integer 
p and q. We will refer to the corresponding Riemann surface as Mp,q- 

Physically, what is happening is that we start probing the target space at large positive 
X, where the LiouviUe field is weakly coupled and the classical target space is a good 
description. Next, we bring the FZZT brane probe into the strong-coupling region x of 
order one, and we find a branch point at x = —1. This branch point is a sign that 
the target space is modified due to strong-coupling effects on the worldsheet. It suggests 
that we analytically continue into the complex x plane, through the branch cut ending at 
X = —1, and into the other sheets of the Riemann surface. Correspondingly, the moduli 



^ This intuition can be made more precise using the equivalent description in terms of the 



Backlund field (p. Here the FZZT brane corresponds to a Dirichlet boundary condition on (f> [23|. 
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space of Y7/LT branes must be enlarged from the complex x plane to the Riemann surface 
■M.p^q- In this way, we obtain M.p^q for the semiclassical target space of the minimal string.il 
A special situation occurs when J^p,q has genus zero. This happens for backgrounds 
without ZZ branes []23|j2^] . Since here M.p,q has genus zero, it can be uniformized by a 



single complex parameter z; i.e. there is a one to one map between the complex z plane 
and M.p,q- Then we can express the values of x and y as polynomials of degrees p and q 
in z 

X = Xy{z) , y = yq{z) (1.5) 

The parameters in these polynomials depend on the closed string background which is 
labelled by the coefficients of closed string operators t = (ti, ^2, • • •)• Of particular impor- 
tance among these parameters are the coefficient of the lowest dimension operator 

r = h (1.6) 

and the worldsheet cosmological constant ^ (in unitary worldsheet theories t = n, but in 
general they are different) . The uniformizing parameter z has the worldsheet interpretation 
as the one point function on the disk of the lowest dimension operator 

z = dr^ = dr(^j ydx^ (1.7) 

Below we will relate this expression to various points of view of the minimal string, in 
particular the connection to integrable hierarchies. 

Most of the discussion in this paper will concern the class of backgrounds where M.p,q 
has genus zero. We expect our main results to apply to the more generic backgrounds with 
arbitrary Riemann surface and we will comment on this below. 

A further useful specialization of the background is to the conformal backgrounds of 
m. These are backgrounds without ZZ branes in which all the closed string couplings t. 



with the exception of the worldsheet cosmological constant have been set to zero. In 
this class of theories explicit worldsheet calculations can be performed, leading to checks 



^ It may seem strange that we started with a one dimensional target space consisting of 0, and 
we ended up with a two dimensional target space consisting of Mp,q- However, the fact that all 
physical quantities depend holomorphically on M.p,q suggests that there is still a sense in which 
the target space is one dimensional. 
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of the results in more general backgrounds. Here the more general expressions (|1.4p , (|1.5|) 
become H (see also 



F{x,y)=Tp{y/C)-T,{x)=0 (1.8) 

and 

x = Tp{z), y = CT,{z) (1.9) 

Here C is a normalization factor which we will determine below, T„(cos^) = cos(n^) is a 
Chebyshev polynomial of the first kind, and for simplicity we have set /i = 1. 

It is natural to expect that this picture of the classical target space is modified only 
slightly when perturbative effects in the string coupling Qs = h are taken into account. 
As we will see, however, nonperturbative effects have important consequences. To study 
the quantum target space, we turn in the first part of section 2 to the nonperturbative 
description of FZZT branes afforded by the dual matrix model. In the matrix model, 
FZZT branes are described by insertions of the exponentiated macroscopic loop operator 
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<i>{x) ~ e^-^^^sC^-A^) = det(a; - M) (1.10) 



into the matrix integral. We use the results of ||3^ to compute the correlator of any number 



of FZZT branes. Taking the continuum limit, we show that the correlators become 

Aid, 




'{lli^ix,,t) (1.11) 



A{x) . ^ 
where 

'il;{x,t) = {^{x)) (1.12) 

is the FZZT partition function, which depends on the closed-string couplings t = ti, t2, ■ ■ ■; 
A denotes the Vandermonde determinant; and dj is shorthand for the action of d = hdj- 
on "tp^Xj.t). (Note that it does not refer to differentiation with respect to Xj.) We argue 
that the denominator can be removed by thinking of the FZZT branes as wavefunctions 
(half-densities), with the result that the branes become fermionic. 

In the classical /i — limit, we show that the FZZT correlators reduce to 




where 



is the semiclassical approximation to the FZZT partition function. We also provide a 
worldsheet interpretation for the various factors in (|1.13|) - the first factor comes from 



annulus diagrams between different FZZT branes, while the semiclassical wavefunctions 
^cK^i) come from the disk amplitude and the annulus between the same brane. More 
generally, we interpret these expressions as a change in the measure of the D-branes \E'(x) 
to a fermion on the Riemann surface Aip,q- 

Using the fact that il^{x, t) is a Baker- Akhiezer function of the KP hierarchy, which is 
actually an entire function of x [|35| , |3^ , it follows that the exact FZZT correlators ( |1 . 1 IJ ) 



are all entire functions of x. This is in spite of the fact that the classical correlators ( |1.13| ) 
are clearly functions on A4p,q- Evidently, the quantum target space differs significantly 
from the classical target space. Whereas the latter comprised the Riemann surface Aip^q, 
the former consists of only the complex x plane! 

In section 2.4, we analyze additional FZZT observables (the quantum resolvents) and 
show that they are also entire functions of x. Finally, we rederive the WKB approximation 
( p..l4| ) using the fact that iIj{x, t) is a Baker- Akhiezer function. We note that the asymp- 
totics exhibit "level crossing" behavior at large negative x. Here, by level crossing we mean 
simply that there is a branch cut along the negative real axis with different values ofilj{x, t) 
above and below the cut. Below, when this approximation of ij{x, t) will be associated 
with saddle points in an integral, we will see that two saddle points exchange dominance 
there. 

In section 3, we illustrate our general arguments with the simplest example of minimal 
string theory, namely the topological (p, q) = (2, 1) model. This is dual to the Gaussian 
matrix model, and we show that the FZZT partition function is expressed in terms of the 
Airy function. By representing the insertion of a D-brane as a Grassmann integral in the 
matrix model, we give a direct and simple proof of the equivalence between the n x n 
Kontsevich model and the double-scaled (2, 1) model with n FZZT branes. That is, we 
show that in the continuum double scaling limit 

with 5' and X nxn Hermitian matrices and r, which can be absorbed in X, is the coupling 
constant of the theory. The eigenvalues of X after an appropriate shift and 




5 



rescaling in the double scaling limit (see below). Using this approach, we see very directly 
how the matrix S of the nxn Kontsevich model is the effective degree of freedom describing 



open strings stretched between n FZZT branes [24 



We continue our study of the (2, 1) model in section 4, focusing now on the effective 



theory on the FZZT brane, which is described by the Airy integral (|1.15|) with n = 1. 
An analysis of this integral using the stationary phase method reveals several new facts. 
To begin, we show how the other sheets of the classical moduli space can be viewed as 
saddle points in the integral describing the FZZT partition function. Therefore, they can 
be thought of as instantons in the effective theory on the brane. We expect this conclusion 
to hold for all values of (p, q). 

A more careful stationary phase analysis of the Airy integral illustrates the general 
mechanism by which the target space is modified nonperturbatively. Exponentially small 
quantities - neglected in perturbative string theory - can become large upon analytic con- 
tinuation, and these large corrections "erase" the branch cuts and monodromies of the 
Riemann surface in the exact answer. The essence of the replacement of the semi classi- 
cal target space M.p,q by the humble complex x-plane is thus what is known as Stokes' 
phenomenon. Generally speaking. Stokes' phenomenon is the fundamental fact that the 
analytic continuation of an asymptotic expansion can differ from the asymptotic expansion 
of an analytic continuation. In our case, the Riemann surface is extracted by working in 
the classical approximation to string theory (thus taking the leading term in an asymptotic 
expansion in Qs = h) and then considering the analytic continuation. Thus, the Riemann 
surface arises from the analytic continuation of the asymptotics. Thanks to the matrix 
model we can study the analytic continuation of the exact nonperturbative answers for am- 
plitudes directly. The fact that the FZZT amplitudes are entire shows that the Riemann 
surface "disappears" nonperturbatively. 

In terms of the saddle point analysis of the field theory living on the FZZT branes. 
Stokes' phenomenon is exhibited in two ways. The first, more trivial way, occurs when 
the parameter x is varied across what is known as an anti-Stokes' line. This can be 
thought of as a first-order phase transition where two contributing saddle points exchange 
dominance. The story is incomplete, however, if we simply consider only the anti-Stokes' 
lines. In addition, there is also a more subtle phenomenon happening along what are called 
Stokes' lines. As x is varied across a Stokes' line, a subdominant saddle abruptly ceases 
to contribute to the exact answer. This phenomenon is most dramatic when we continue 
to vary x and the missing saddle becomes the dominant saddle, even though it is still not 
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contributing to the integral. In terms of the path integral describing the effective theory on 
the brane, what is happening is that one simply cannot deform the contour of integration 
to pass through that saddle. We will discuss this in more detail in Appendix B. 

In section 5, we comment on the issues involved in generalizing to other backgrounds. 
Among other things, we show using the FZZT partition function that not all values of 
(p, q) correspond to nonperturbatively consistent backgrounds with a double-scaled matrix 
model that is bounded from below. We deduce a bound 

sin— >0 (1.16) 
P 

that must be satisfied in order for the corresponding background to exist. For instance, 
when p = 2 only the [p, q) = (2, 2m — 1) models with m odd exist nonperturbatively, while 
the models based on the unitary discrete series with q = p+1 never exist nonperturbatively. 

Finally, section 6 contains a possible analogy with the work of on the physics 
behind black hole horizons and possible implications for the topological string approach of 
3^ . In appendix A, we review the Lax formalism of minimal string theory (the operators 



P and Q and the string equation [P, Q] = h), and we present new results concerning the 
geometrical interpretation of its classical limit. Appendix B, as mentioned above, contains 
a brief review of Stokes' phenomenon along the lines of [^, while appendix C contains 
the results of a numerical analysis of (p, q) = (2, 5). 



2. The Quantum Target Space: FZZT Branes in the Matrix Model 

2.1. FZZT correlators at finite N 

As is well-known, (p, q) minimal string theory possesses a dual matrix model descrip- 
tion. For p = 2, the dual matrix model consists of an x Hermitian matrix M with 
potential V{M) and coupling g, 

Z{g) = j c^Me-iT'^^(^) (2.1) 

while for p > 2 one needs to use an analogously defined two-matrix model (for recent 
discussion of the two- matrix model and references, see ||4C1|| ): 

Zig) = I dMdM .-^{^^viMHTrwm-T.MM) ^^.2) 
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Here the measures dM and dM include a factor of the volume of U{N). 

In the matrix model, macroscopic loops are created by insertions of the operator 

W{x) = -^Tr \og{x - M) (2.3) 

in the matrix integrahi For instance, the large N limit of (yV{x)) corresponds to the FZZT 
disk amplitude $ (up to a polynomial in x), and the matrix model resolvent 

R{x) = d^{W{x)) = 1 (^^^) (2.4) 

corresponds to y{x) (again up to a polynomial in x). 

The full FZZT brane obviously does not correspond to a single macroscopic loop 
in the worldsheet. Rather, we must include contributions from worldsheets with any 
number of boundaries. This is accomplished by exponentiating W{x), whereby the full, 
nonperturbative FZZT brane is represented by a determinant operator 

gW(x) ^ ^g^^^ _ ^) ^2.5) 

in the matrix model. We can also write this determinant as a Grassmann integral over 
complex fermions Xi 

det{x-M) = [ dxdx^ e^^^^-^)^ (2.6) 



In [p!8| - pT| , p6| , the matrix M of the one- matrix model was interpreted as describing the 
(bosonic) open strings stretched between the condensed ZZ branes in the Fermi sea. 
Meanwhile, the Xi are taken to represent fermionic open strings stretched between the 



FZZT brane and the ZZ branes 26 



Now consider the correlation function of any number of FZZT branes, which non- 
perturbatively is given by a product of determinants. Amazingly, this can be explicitly 
evaluated in both the one and two matrix models . The answer is 

ndet(..-M)\ = !5!i(^gt^ (2.7) 



^ In the two-matrix model, there is another loop made out of M. It corresponds to the "dual" 
FZZT brane and classically is related to the loop ( ^.3|) by a Legendre-type transform ||40| ]. We 
will discuss the interpretation of this dual loop in section 6. 
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Here A(x) = ni<,( 

the Vandermonde determinant, Pk{x) are the orthogonal 
polynomials of the matrix model (or bi-orthogonal polynomials associated to M in the 
two-matrix model) with leading coefficient 1, and the indices i and j in (|2.7D run between 



1 and n. The simplest case of the general formula (|2 . 7|) is the FZZT partition function. 
This is given by a single orthogonal polynomial: 

{dei{x - M)) = Pn{x) (2.8) 



Below, we will take the continuum limits of (|2.7|) and ( p.8|) , and we will see how the 
perturbative loop correlators can be recovered. 

Before we proceed, let us briefly mention an interpretation of the FZZT correlator 
( p.7|) that will be useful in the next subsection. First, we need to write the LHS of ( |2.7| ) 
more compactly in the following way 

JJ det(xi -M) \ = (det(X ®lN-In®M)) (2.9) 
u=l / 

where In and denote the N x N and n x n identity matrices respectively, and X is 
understood to be an n x n Hermitian matrix with eigenvalues xi, . . . , Xn- Now notice that 
if we square det{X ® In — In ® M), multiply by e~'^^^^^^^^~^^^^'^^\ and integrate over 
X and M, we obtain the (A^ + n) x (N + n) matrix integral with no insertions of FZZT 
branes, i.e. 

j dXdM e--^^'''^^^+^'''^''^^det{X ® In - In® Mf = j dM e--^^^''^^'^ (2.10) 

where M is an (A^ + n) x (A^ + n) Hermitian matrix. (As before, the integration measures 
include factors of the volume of the relevant unitary group.) The meaning of ( p. 101 ) is 
that the FZZT creation operator det(X ® In — In ® M) acts as a kind of wavefunction 
(half-density) on the space of Hermitian matrices X. 

The motivation for interpreting det(X ® In — In ® M) as a wavefunction in X-space 
is that it allows us to think of the FZZT branes as fermions. To see this, recall that the 
measure for an integral in X space is 

n 

dX = dU^dx^ /\{xf (2.11) 

i=l 

where [/ is an n x n unitary matrix (and the measure is such that J dU = 1). Hence, a half- 
density {dXY^"^ carries with it a factor of A(a;), which is precisely what is needed to cancel 
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the denominator of (|2.7| ).B (Put differently, the factor A(x) plays a role analogous to that 
of CO cycles in vertex operator algebra theory, by enforcing the correct statistics of the de- 
terminant operator.) This leaves the numerator of ( |2.7| ), which is obviously antisymmetric 
under interchange of the x^'s. Therefore, the FZZT branes become fermionic. 

2.2. FZZT correlators in the continuum limit 



Now let us take the large double-scaling limit of (|2.7| ) to obtain the D-brane cor- 
relators of minimal string theory. For simplicity, we start with the n = 1 case ( |2.8| ). As 
in (|2.1C1| ), to have a well-defined scaling limit we must consider not the determinant, but 
rather the following operator [p5|p| 



^(x) = ^=e-^(^)/2Met(x - M) (2.12) 
VhN 

where V{x) is the matrix model potential, and hjy is a normalizing constant. (Some rigor- 
ous results on the double-scaled limit of the orthonormal wavefunctions have been derived 
in [0.) This converts the orthogonal polynomials in ( p.7|) to orthonormal wavefunctions 
with measure dx. Then the FZZT partition function in the double-scaling limit is given 
by a function of x and the background closed-string couplings t = (ti, t2, ■ ■ 

{^{x))=i(j{x,t) (2.13) 

which is characterized by the requirement that it satisfy the differential equations 

Qip{x,t) = xijj{x,t), Pipix.t) = nd^ijj{x,t) (2.14) 

with P (X d*^ + . . . and Q (x d^ + . . . differential operators in d = hdr^ (Note that the 
derivative d is taken at fixed values of x, tj>i.) P and Q are known as Lax operators, and 
they are determined by the string equation 

[P,Q] = h (2.15) 



^ The annulus diagram is the logarithm of the n = 2 version of ( |2.7| ). In p^ , where this 
diagram was calculated, it was pointed out that the term associated with the denominator of 
(|2.7| ) is independent of the coupling constants and therefore could be removed. Here we see a 
more geometric way of deriving this fact. 

^ It is common in the literature to denote the lowest dimension coupling by x. We denote it 
here by r. 
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In appendix A, we review the properties of P and Q and present new results concerning the 
geometric interpretation of these operators in the classical limit. For a more pedagogical 
introduction to the Lax formalism and integrable hierarchies, see e.g. |T0| , ^ . 

Note that the differential equations ( |2.14| ) do not specify uniquely. In non- 

perturbatively consistent models, we will see below that this ambiguity can be completely 
fixed, in part by the boundary condition that il}{x,t) be real and exponentially decreasing 
as a; — s> +00. 

In the literature on integrable systems, the function ilj{x,t) is referred to as the "Baker- 
Akhiezer function" of the associated KP hierarchy defined by P and Q. Here we see that 
it has a simple, physical interpretation in minimal string theory as the FZZT partition 
function. For our present purposes, the most important property of the Baker-Akhiezer 
function is the non-trivial fact that it (along with all of its derivatives) is an entire function 
of x.i We will see momentarily that this has dramatic consequences for the quantum moduli 
space of FZZT branes. 

In the double-scaling limit, it is not difficult to show that the general FZZT correlator 
( p77| ) becomes 




(2-16) 

i=l 



\i=l 

where the notation dj is shorthand for the action of d = hdr on i(^{xj, t). To derive ( |2.16| ) 
use the fact that the increase of index on Pk becomes a derivative with respect to r to 
leading order in the double-scaling parameter e. Then as e 0, only the Vandermonde 
determinant of derivatives with respect to the index survives. We conclude from ( |2.16| ) 
that the correlator of any number of FZZT branes reduces, in the continuum limit, to a 
product of Baker-Akhiezer functions il^{xi, t) and their derivatives. 

2.3. Comparison with the semiclassical limit 

Having obtained the exact D-brane correlators ( |2.16|) , it is straightforward to take 
their semiclassical limit and show how the perturbative loop amplitudes can be recovered. 



^ This fact was proven in |35,36| by writing the string equation in an equivalent form as an 
equation for a flat holomorphic vector bundle on the space of x, tfe. The connection on this vector 
bundle is polynomial in x. The Baker-Akhiezer function is used to make a covariantly constant 
frame. From the equation (^ — Ax)'^ = 0, where ^ is the frame, it follows, via the path-ordered 
exponential, that ^ is entire in x. 
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For this, we will need a result from subsection 2.5, namely that as — > 0, the Baker- 
Akhiezer function becomes an eigenfunction of d with eigenvalue z. (As discussed in the 
introduction, the global uniformizing parameter z exists only in the backgrounds without 
ZZ branes.) Therefore, the semi classical limit of ( |2.16| ) is simply 




The worldsheet description of the various terms appearing in (|2.17|) is as follows. Recall 
that we can think of the FZZT creation operator \E'(x) ~ e^(-)/^ as the exponentiated 
macroscopic loop operator. Then the first factor is the exponentiated contribution 
of the annulus diagrams with the ends of the annulus ending on different branes. This is 
consistent with the explicit worldsheet calculation in conformal backgrounds which leads 
to the connected annulus amplitude p6| 

(T^(x)T^(x'))c,ann = log (2-18) 

The other diagrams that contribute at this order in h are the disk diagram (|1.3| ) and the 
annulus diagram with the two ends on the same braneB 

hm ^(Vr(x)Vr(x'))c,ann = hm ^ log { \~ \ = -^logdM^) (2-19) 
These diagrams combine to give the WKB wave functions in ( |2.17| ): 

/x{z) 
ydx (2.20) 

1 



\f%x[z) 



As is common in WKB wavefunctions, the prefactor j{z^ is a one loop correction. In our 
case it arises from an open string loop which is the annulus diagram. Finally, note that 
higher genus diagrams are suppressed in the /i ^ limit in (|2.17D 



^ The factor of a half comes from the fact that the open strings are ending on the same brane. 

^ We would like to stress that the leading order expressions ( |2.17D , ( p.20| ) are correct in any 
background without ZZ branes and not only in the conformal backgrounds. The only fact that 
is needed is that the Riemann surface A^p,<j can be uniformized by the complex parameter z; i.e. 
that it has genus zero. We will return to this WKB wavefunction in section 2.5. 
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The difference between the classical result ( p.l7|) and the exact quantum result ( |2.16| ) 
is at the heart of our analysis. The classical answer is obviously defined on a multiple 
cover of the complex x plane, since for the same x, there can be p diff'erent values of z. 
On the other hand, since il){x^t) and its r derivatives are entire functions of x, the exact 
expressions (|2.16D for the FZZT correlators are actually entire in the complex x plane. 
In other words, there are no branch cuts or other singularities that necessitate analytic 
continuation to other sheets. Apparently, the semiclassical target space M.p,q disappears 
when one takes nonperturbative effects into account! 

In the next subsection, we will check our picture of the quantum target space by com- 
puting the quantum resolvent and showing that it is also an entire function of x. However, 
we would like to first mention another perspective on the semiclassical correlator ( |2.17| ) 
and how this is modified in the exact answer. Recall that the FZZT brane could be thought 
of as a half density multiplied by (dx)^/^. Thus the annulus factor f{z) can be interpreted 
as a measure factor implementing a transformation from x to z. By the same token, we 
can think of the correlator of n FZZT branes as a half-density multiplied by ((iX)^/^ where 
X is an n X n matrix with eigenvalues Xi (see ( ^.11| )). Then the transformation of this 
half-density to Z-space, with Z an n x n matrix whose eigenvalues are Zi, must include a 
factor of the Jacobian 

dZ 



dX 



\{f{^^) (2.21) 

i=l 

But according to the discussion above, this is precisely the contribution of the annulus to 
the correlator! Thus we have shown that 

hm (^nxi)^ {dXf^ = e^'^^^^HdZy/^ (2.22) 

with $(z) as in ( p.20|) . Evidently, the classical correlators reduce to extremely simple 
expressions in 2;-space. 

These simple expressions suggest there should be an equally simple formalism under- 
lying the classical theory. One possibility was alluded to above, namely that instead of 
thinking of these correlators as half densities we can equivalently think of them as fermions. 
Then (|2.22| ) indicates that in the semi- classical limit, there is a sense in which the FZZT 
branes are actually fermions on the Riemann surface M.p,q- (Remember that the z plane 
covers M.p^q exactly once.) Such fermions are common in the matrix model literature (for 
reviews, see e.g. [^ ^ , ^ , ^ ) . However, our general discussion suggests that this simple 
picture cannot be correct in the full nonperturbative theory, in which the Riemann surface 
M.p,q is replaced by the complex x plane. We return to this point at the end of section 6. 
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2.4- The analytic structure of the quantum resolvent 

Although we have shown that physical observables involving the determinant operator 
are entire functions of x, it remains to be seen whether the same is true for the resolvent 
R{x) defined in ( |2.4| ). To all orders in h the resolvent is expected to exhibit monodromy 
and have various branch cuts in the complex x plane. But in the dual string theory, the 
resolvent (or rather its integral) corresponds to the vacuum amplitude of a worldsheet 
with one boundary and an arbitrary number of handles. Thus we might expect that 
nonperturbative effects drastically modify the classical resolvent, just as they modified the 
classical determinant correlators (|2.17| ). 

For simplicity, we will limit ourselves to the one-matrix model, which describes the 
theories with p = 2. Then the matrix integral defining the resolvent can be easily reduced, 
using the determinant formula ( |2.7| ), to a single integral in terms of orthogonal polynomials: 

R{x) = r ^ (2-23) 



oo 




i;M-i{X)i^NiX) - i^N-iWN{X)) (2.24) 



where 



with il)k{X) = —^e ^^^^^"^^ Pk{X) the orthonormal wavefunctions of the matrix model and 
tp' = -^''P- As we discussed above, tpNiX) becomes the Baker- Akhiezer function 'tjj{Xjt) in 
the double-scaling limit. Then the exact, double-scaled resolvent iM 

R{x) = r (2.25) 

J_oo X-X 

with 

pn{X) = Ah^ {dMX, t)c>AV(A, t) - ^(A, t)a,c>A^(A, t)) (2.26) 

where A denotes some overall numerical factor which will be irrelevant for our purposes, 
and r corresponds to the lowest-dimension coupling as below ( |2.14| ). From the expression 
for the double-scaled resolvent, it is clear that we can think of Ph{X) as defining the 
quantum eigenvalue density. 



^ As is usual when defining the continuum resolvent, one might have to impose a cutoff on the 
integral at —A; this does not affect the arguments below. 
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Since the Baker- Akhiezer function is an entire function of A, the resolvent will be 
everywhere analytic, except along the real axis, where it suffers from a discontinuity 



R{x + ie) - R{x -ie) =2nipn{x), xeR (2.27) 

Contrast this with the classical resolvent, which is discontinuous only along a semi-infinite 
cut. The discontinuity (|2.27| ) suggests that we define two resolvents, R+{x) and R-{x), 
which are obtained by analytically continuing R{x) through the real axis from either the 
upper half plane or the lower half plane, respectively. Explicitly, we define 

R±{x) = / ^ (2.28) 

JC± X — A 

where the contour C_|_ (C_) travels below (above) x and satisfies the same boundary 
conditions at infinity as the original contour in ( p.25|) . Then we have 

o/^N f^+(^) forlmx>0 , . 

^^^^ ~ \ i?_(x) for Imx < ^ 

and also 

R+{x) - R-{x) = 2mpn{x) (2.30) 

for all a; G C. Given the definition ( |2.28| ), it is clear that both R±{x) are entire functions 
of X. 

Finally, let us consider the classical limit h 0. In this limit, the resolvent must 
reduce to the classical resolvent, which solves the factorized loop equation and has a semi- 
infinite branch cut along the real x axis. Therefore, according to ( |2.29D , R+{x) has the 



correct classical limit for Imx > 0, while R-{x) has the correct classical limit for Ima; < 0. 
Analytically continuing the classical limits of either R^{x) or i2_(x), we find the branch 
cut and the second sheet of the Riemann surface. Note that it is essential first to take the 
classical limit (i.e. drop the nonperturbative corrections), and only then to analytically 
continue the resolvent. Otherwise, we will not find the second sheet, since R±{x) are both 
entire functions of x. 

To summarize, we have seen that it is impossible to define globally the quantum 
resolvent R{x), due to the discontinuity on the real x axis. Instead, we can define through 
analytic continuation two resolvents R±{x), both of which are entire in the complex plane. 
So for the resolvent, just as for the determinant, the Riemann surface disappears at h ^ 
and is replaced with the complex plane. To recover the Riemann surface, we must first 
take the classical limit of the resolvents, and then analytically continue. 
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2. 5. More on the FZZT partition function 

We have seen above how in the double-scaUng hmit, all of the observables involving 
the FZZT brane reduce to products, derivatives, and integrals of a single quantity, the 
FZZT partition function i(j{x,t). Thus it makes sense to study this object in more detail. 
Although we do not have a general formula for il)[x^t) (see below however, where we 
study the example of the Gaussian matrix model), we extracted its asymptotic behavior 
at small h in (|2.20| ) using worldsheet methods. Here we would like to rederive the WKB 
approximation 

^ ^d{z, t) = {d,x{z, t))-^/h-^T' (2.31) 

starting from a completely different point of view, namely the fact that the FZZT partition 
function is a Baker- Akhiezer function of the KP hierarchy. We must demonstrate that 
( p.31|) satisfies ( p.l4|) in the h limit. This was first shown in for p = 2 (together 
with the interpretation in terms of Riemann surfaces). We now give a simpler, but equally 
rigorous, proof of this result. 

The first step in the proof is to act on with d — hd^ (at fixed x). To leading order 
in h, this gives 

= / d^yldx + 0{h) (2.32) 

^Cl JxQ 

We can simplify this by writing y = y{x{z,t),t) and converting the derivative at fixed x 
to one at fixed z: 



/ dry\^dx = / (dry\^~ dxy\^drx\]dx 

Jxq Jxq 

/ [dTy\^d;,x\^- d^y\^drx\)dz 



(2.33) 



where we have used dzy\^ = dxy\^dzx\_^ in the second equation. We recognize the integrand 
in the second equation to be the Poisson bracket of x and y. Using the freedom to shift 
zq ^ 0, together with the fact that x and y are given by 

x{z,r)=Q{d = z,r)\^^^, y{z,T) = P{d = z,r)\^^^ (2.34) 

and must therefore satisfy the genus zero string equation 



{x, y} = drxd^y - dryd^x = 1 
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(2.35) 



(see appendix A for a proof of this), we conclude that 

drvl dx = z (2.36) 



Using ( p.36| ) in ( |2.32| ) we readily see that the classical Baker-Akhiezer function is an 
eigenfunction of d with eigenvalue z in the classical /i — > limit. Note that, as mentioned 
in the introduction, this identity can also be proven using worldsheet techniques for the 
special case of the conformal background [^. The advantage of the derivation we have 
given here is that it is valid in every background where the uniformizing parameter z exists 
(i.e. backgrounds without ZZ branes, in which case M.p,q has genus zero). 

At the next order in h we need to consider several terms. The first step is to expand 
the operators Q = Qo + hQi + ■ ■ ■, and P = Pq + hPi + ■ ■ ■, where all derivatives are on 
the right hand side. In appendix A, we show that 

Qi = ^drd,x\^^^, Pi = ldrd,y\^^^ (2.37) 



where x{z,t) and y{z,T) are the classical expressions ( p.34| ). (Except when specified ex- 



plicitly, the derivatives of x are taken when x{z, r) is considered as a function of z and r. 
Thus dr here is taken at fixed z,tjyi.) When d acts on the exponent of the wavefunction 
( p.31|) it gives back z as in ( |2.36|) . In addition, we need to consider two more terms at this 
order in h. The first contribution arises from a second derivative of the exponent. This 
gives a term of the form 

^afx {drz\,) (2.38) 

where the factor of ^d'^x comes from selecting the two derivatives in Qo which are acting 
twice on the exponent of the wavefunction, and then evaluating the rest of the derivatives 
using the classical result. The second term appears when the derivatives of Q act on the 
prefactor of ( 2.31| ). This leads to a term of the form 



~d,x dAlogd,x\r)\, (2.39) 

where again the factor dzX selects the derivative in Qo that is acting on the prefactor of 
( p.31|) . So finally we obtain that 



(Q - x)<f,i = h ( -dix d^z\, - -d,x dr (loga,x|,)|^ + Q^ ] + 0{h') (2.40) 
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We now use Ot-zI = — to simplify these terms. In particular, we have 

dr (loga^x = ^ = ' 2.41 



Then using (|2.41[) and ( p.37| ) we find that all terms of order h in ( |2.40| ) cancel. 



Computing the action of P on takes a little more work. Using again the results 
above, one can show that 

(P - hd.)^., = h ( ^^^-y^^"" - - ^-^^-^-^ + dMdzydrX-l)\ ^ ^(^2) 

2 V 9zX [OzxY J 

(2.42) 

By applying the genus zero string equation ( |2.35|) and its derivative with respect to 2, we 
see that the terms in parentheses all cancel, confirming that 

p^d = riB^'^ci + o{n^) (2.43) 

This completes our proof that '^ci is indeed the leading-order WKB approximation to the 
Baker- Akhiezer function. 

Let us also offer the following non-trivial consistency check of the semiclassical approx- 



imation (|2.31|) . This approximation clearly suffers from a p-fold ambiguity, corresponding 
to which branch of y{x) and which value of z in x'{z) we choose. The correct branch 
is chosen at large positive x by demanding that y{x) be given by its physical sheet as 
X — i> -|-oo. For large |x| in the first sheet we have 

y^2p~^Cxp (2.44) 

where the real constant C was determined in Its sign is 



r] = sign(C) = —sign (sin {qTt/p)) (2-45) 

This means that up to a power of x 

'(/'(x, t) ~ exp ( ryCx p j (Imx = 0, Rex^+oo) (2-46) 

with C real and positive. We expect this semiclassical approximation to be valid every- 
where at large |x|, except on the cut along the negative x axis. But then the fact that 
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il){x^t) is entire means that as we cross the cut, the asymptotic behavior of i(j{x,t) must 
change from 

/ ~ p+q \ 

ip{x^t) ~ exp (rjCx p j (Imx > 0, Rex — > — oo) (2.47) 
above the negative real axis, to 

tjj{x,t) exp (^rjCe~'^''''(^'>x^^ (Imx < 0, Rex ^ -oo) (2.48) 

below the negative real axis. Slightly above and below the cut, both contributions ( |2.47D - 
( p.48|) are present. In order for this to be consistent with the semiclassical approximation, 
the first contribution must dominate above the cut, while the second must dominate below 
the cut. Fortunately, this is guaranteed by the sign of 77 = sign(C) ( |2.45[ ). It is very 



satisfying to see how the semiclassical approximation, the level crossing behavior, and the 
worldsheet calculation of C all fit together so consistently. 

Finally, it is worth mentioning that this level crossing behavior is an example of 
Stokes' phenomenon. We review Stokes' phenomenon in appendix B, and we will discuss 
its implications for the quantum target space in much greater detail in section 4. 



3. An Example: The Simplest Minimal String Theory and its FZZT Brane 

3.1. FZZT correlators and the Kontsevich model 

Here we will illustrate the ideas of the previous sections using the example of the 
(p, q) = (2,1) model, also known as topological gravity. Since this theory is dual to the 
Gaussian matrix model, it allows us to make quite explicit some of the general formulas 
derived above. Along the way we will encounter a new point of view on the relationship 
between the Kontsevich matrix model and the double-scaling limit of matrix integrals. 

The (2, 1) model is represented in the matrix model by the integral 

Z{g)^ j dMe--^^'^" (3.1) 



Another argument that the Baker- Akhiezer function must exhibit Stokes' phenomenon, based 
on the behavior of eigenvalue distributions, was given in sec. 3.8 of [p5|. This argument is related 



to work of F. David on nonperturbative stability [44| and is also in accord with the discussion of 



nonperturbative stability at the end of section 5 below. 
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with M an N X N hermitian matrix. An FZZT D-brane insertion is represented by the 
integral 

(det(a; - M)) = j dM det(a; - M)e-^^'^^ (3.2) 

Using (|2.6|) , we can write this as an integral over the matrix M and N fermions Xi- Then 
we can easily perform the Gaussian integral over M in to find the effective theory of the 
fermions 

(det(.-M)>^/drfe-fx'x>'--'^ (3.3) 

Note that we started in with A^^ grees of freedom, the entries of M. After gauge 
fixing they are reduced to the eigenvalues of M, whose number is A^. Now we have order 
A^ fermions, but their effective theory - which is still invariant under U{N) - depends only 
on a single variable x^X- To make it more explicit we replace (|3.3|) with 



(det(a;-M)) = j dsdxdx^ e-i-'+(-+-)x^x _ ^ 1_ ^ ^ i^)^ (3.4) 



and view s as an effective degree of freedom on the FZZT brane. The final expression as 
an integral over s is similar to the starting point ( |3.2| ). The matrix M is replaced by a 
single variable s and the determinant is replaced with [x + is)^ . 

We recognize the RHS of ( |3.4| ) as the integral representation of the Hermite polyno- 
mials: 

(det(x-M)) = (|)*if^(x^) (3.5) 

Since these are the orthogonal polynomials of the Gaussian matrix model, this confirms 
explicitly in this example the general result ( |2.8|) . 

It is trivial to generalize this discussion to n FZZT branes. The partition function of 
n FZZT branes is given by 

(det(X ® - ® M)) = -i^ / dMdxdx^ ^-^TrM-+xHxm^-i.m4)x ^3_g) 

with X an n X n matrix and Xaj, Xaj fermions transforming in the bifundamental repre- 
sentation of U{n) X U{N). Integrating out M and integrating back in an n x n matrix 5", 
we find (after dropping an overall factor) 

(det(A: ^lN-In^M)) = J dS det{X + iS)^e~^^'^^ (3.7) 
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In the large limit with g ~ the eigenvalues of M become localized in the 

interval (— "\/2, "\/2) along the real axis. The double-scaling limit then corresponds to 
zooming in on the end of the eigenvalue distribution, while simultaneously bringing the 
two saddles of ( |3.4| ) together. For example, for n = 1 the double-scaling limit of the FZZT 
partition function (|3.4| ) is 

x^y/2(^l + ^e^h-^/^x^ , Ng^l- e^h'^/^T, s ^ (z - es) , N ^ e'^ (3.8) 



with e — s> 0. Here r is the lowest-dimension coupling in the continuum theory. Then 
becomes (after dropping overall numerical factors) 

/oo 
^i,,3+,n~^/3^,+r)s (3.9) 
-oo 

We recognize this as the Airy integral; therefore the FZZT partition function is simply 

ij{x,t) = Ai{{x + T)h-^/^) (3.10) 

There are a few things to note about this result. 
1. The FZZT partition function ( |3.1C1| ) clearly satisfies 

Qtp = xtp, Ftp = hd^tp (3.11) 

with Q and P given by 

Q = (f + r, P = = d (3.12) 

(These operators obviously satisfy the string equation [P, Q] = h.) This confirms, in 
this example, that the FZZT partition function is the Baker- Akhiezer function of the 
KP hierarchy. 



2. The Airy function (|3.10| ) is an entire function in the complex x plane. On the real axis, 
it is oscillatory for x < — r (where classically there is a cut) and decays exponentially 
for X > —T. Therefore, although there appear to be two FZZT branes with the same 
X semiclassically (corresponding to the different sheets of ^^2,1)7 we see that the fully 
nonperturbative FZZT branes depend only on x. The Riemann surface disappears 
nonperturbatively and is replaced with only its physical sheet. 
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3. There are, of course, two linearly independent solutions to the equations The 



other solution is the Airy function Bi. We see that it does not correspond to the phys- 
ical FZZT partition function. Indeed, this solution behaves badly in the semiclassical 



regime x — > +00, where it grows exponentially. In terms of (|3.9|) Bi corresponds to 
another integration contour. 



Now consider the analogous double-scaling limit for the general FZZT correlator ( |3.7|) . 
In this limit, we find 

which is, of course, the n x n Kontsevich model (for a review of the Kontsevich model 
and topological gravity, see e.g. [0). Through the use of the fermions, we have obtained 
a rather direct route from the Gaussian matrix model to the Kontsevich model. We also 
see quite explicitly how the matrix 5" of the Kontsevich model is the effective degree of 
freedom describing open strings stretched between n FZZT branes, an insight obtained in 

Note that we can also think of the FZZT correlator ( p.l3|) as a perturbation around 
the closed-string background corresponding to {p,q) = (2, 1). This is a trivial statement 
at finite A^: it simply means that the insertion of determinants at positions xi, . . . ,Xn is 
equivalent to a certain shift in the matrix model potential. This shift can be obtained 
by writing det(a;i — M) as e"'"''^"^^^*"-^^) and expanding the logarithm at large Xi. In 
the continuum limit, essentially the same story holds, except that now we must expand 
log(x — M) in the basis of scaling potentials Wk{M) for a single cut model. For a cut 
between —V2 and -\/2, these potentials take the form H 



r- N 1/2' 



(3.14) 



-I + 



where [ ]+ indicates that we expand in powers of 1/M and keep only positive powers of M. 
Thus as A ^ 00, each determinant insertion det(x — M) can be viewed as a modification 
of the potential 

00 

Tr V'{M) ^ Tr V'{M) + ^tafc+iTr M/^(M) (3.15) 

fc=0 
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where the couphngs are given by 

9-(A;-l)/2 / 9^/9 \~^/^ 



This formula for t2k+i can be verified by, e.g. writing the scahng potentials ( p. 14]) as 
contour integrals around infinity and then performing the sum ( p.l5| ) explicitly. 

In the double scaling limit we zoom in on x ~ a/2 as in ( |3.9|) (we set h = 1 and r = 
for simplicity). Summing over z = 1, . . . , n, we find that the t^'s reduce to 



- 2fc + _ 

1=1 



'-jt-^'-'" (3.17) 



Therefore the (2, 1) model with n FZZT branes can be thought of as the closed-string 
background with the couplings (|3.17| ) turned on. 



To identify properly the precise value of the closed-string partition function, we must 
also take into account the fact that ^closed (i) ^ 1 as t ^ (which is the same as Xi 00). 
In this limit, with Xi +00, the double-scaled FZZT correlator reduces to the WKB 
approximation of the matrix Airy integral ( p.l3|) . Thus we must divide by this quantity 
to extract the closed-string partition function. In other words, we have shown that in the 
double-scaling limit, 

g-TrXV23^^et(X ®lN-In® M)) C(X)Zelosed(t) (3.18) 

with 

C{X) = e-2TrX^/V3 j e-T>-v^^^ (3.19) 
Equating (|3.13| ) and (|3.18| ), (and setting /I = 1, r = 0) we arrive at the relation 



J dS eTr(iSV3+iXS+2X3/2/3) 
^closed (t) = J^g^-TrVXS^ ^^'^^^ 

Note that by shifting the 5" integral in the numerator, we can rewrite this as 

2.closed(^) - J^^g-TrZS2 (3-21) 

with 

Z = VX (3.22) 



23 



Equation ( p.21 ) is the way the relation between the finite n Kontsevich model and topolog- 
ical gravity is usually stated. Here we have rederived this fact directly from double-scaling 
Gaussian matrix model. 

There are a few interesting things to note in our derivation. First, the normalizing 
factor C{X), being the WKB approximation to the (matrix) Airy function, is not an entire 
function of the eigenvalues of X. For instance, it has the simple form (2y^a;)~^/^e~^^^'^^/^ 
when n = 1, and this clearly has branch cuts in the complex x plane.0 This explains why 
the usual relation (|3.21| ) between the closed-string partition function and the Kontsevich 
integral suffers from branch cuts as a function of (the eigenvalues of) X. On the other 
hand, we see that the combination C{X)Zciosed{t), being the matrix Airy integral, is an 
entire function of X , even though the separate factors are not. 

The second point worth mentioning is the interpretation of the quantity Z = \fX (or 



Z = \Or+T for nonzero cosmological constant) that emerged naturally in our derivation. 
This quantity also featured in the work of [^, where it corresponded to the boundary 
cosmological constants of n FZZT branes. In order to compare with the results of ||24|| , 
one needs to keep in mind that in |j2^ the Liouville coupling constant was taken to be 
= 1/2, while here we are assuming 6^ = 2. Thus, their boundary cosmological constant 



is equal to our dual boundary cosmological constant ]1b = ^/fJ-B + t = /J-Blthere- With 
our definitions, the FZZT brane labelled by /xg is the one corresponding to det(/xs — M) 
in the double-scaled Gaussian matrix model. This corresponds to treating the worldsheet 
boundary interaction hb^^'^ as a non-normalizable operator. From the point of view of 
Liouville theory with 6^ = 2 it is more natural to consider FZZT branes as a function of 
JiB, which corresponds to treating the worldsheet boundary interaction JiBe^'^ as a non- 
normalizable operator. The expectation values of the FZZT branes with these two choices 
are related, at the classical level, by a Legendre transform. When we quantize open strings 
on the FZZT brane it looks like we have a choice of which operator to fix. These two 
choices amount to different quantization prescriptions for the open strings, analogous to 
the different choices in other AdS/CFT examples It seems that the open string field 
theory of corresponds to considering fluctuations of JIb- This ends up performing a 
Fourier transform between the result for fixed JIb, which is a simple exponential, and the 
result at fixed /U^, which is given by the Airy integral. 



Note that this is consistent with the asymptotic expansion at large x of our general result 

dm). 
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3.2. The quantum resolvent in the Gaussian matrix model 

The Gaussian matrix model also provides a good setting for the discussion of the 
quantum resolvent in section 2.4. Substituting ( p.lO|) into (|2.26|) , we find the quantum 
eigenvalue density 



(A) = h^/^Ai'{{X + T)h-^/^f - h-^/^iX + T)Ai{{X + T)h-^/^f (3.23) 



One can check that this is everywhere positive on the real axis. Using the asymptotic 
expansions 

f ^ 1/4 e~^/^^^^^ |arg(x)|<7r 
we see that the classical limit of the eigenvalue density on the real axis is as expected: 

limp,(A)= ^^/^ (3.25) 



A<-r 



Now consider the quantum resolvent. Since pfi is positive on the real axis, the quantum 
resolvent ( |2.25|) is indeed discontinuous across the entire real axis. Combining (|2.3C1| ) with 
( p.24|) , we see that the resolvents and R- only differ by a small, nonperturbative 
amount in a wedge around the positive real axis (we now set r = for simplicity). 

R+ix) - R-ix) ~ — e"^"'^', I arg(x)| < - (3.26) 
Ax 3 

Let us call this wedge region I. By the same argument, the resolvents differ by a large 
amount in the wedge ^ < | arg(x)| < tt, which we will call region II. In other words, the 
small nonperturbative quantity ( |3.26| ) in region I becomes a large nonperturbative quantity 



in region II. In region II, the resolvent Rj^ has the correct classical limit in the upper half 
plane, while the resolvent R- has the correct classical limit in the lower half plane. 

To find the second sheet of the Riemann surface, we must first take the classical limit 
of (-R-) in the union of region I and the upper (lower) half plane. Only by dropping the 
nonperturbative corrections does the branch cut at x < appear. Then we can analytically 
continue through this cut to find the second sheet. 
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4. The Effective Theory on the Brane 

We have seen in the previous two sections how the Riemann surface disappears non- 
perturbatively, with the FZZT partition function being an entire function of the complex x 
plane. Here we would like to understand this nonperturbative modification in more detail, 
from the point of view of the effective theory on the FZZT brane. We will limit ourselves 
to the simplest case of (2, 1), in which case the effective theory (|3.13| ) on n FZZT branes 
is the n X n Kontsevich model [^. For simplicity, we will consider only the case n = 1. 



The semiclassical approximation of ^ — corresponds to the saddle-point approxima- 
tion. For the FZZT partition function ( |3.9| ), there are two saddle points in the s integral, 
located at 

(s) = i^-^/^v^ (4.1) 

(For simplicity we set r = 0.) Therefore, there are two distinct branes for each x, semi- 
classically. The moduli space of branes becomes a double cover of the x plane, as we saw 
in the introduction. 

Contrast this now with the quantum theory. Here we must integrate over s; i.e. we 
must study the quantum dynamics of the theory on the brane. The subleading saddles in 
the integral over s are instantons in the theory on the brane. As is always the case with 
instantons, one must sum over them in some prescribed fashion. The result of this process 
is that the exact, nonperturbative FZZT partition function becomes an entire function of 
X. For (2, 1) it is the Airy function ( |3.1C1| ). 

The crucial point is that instead of exhibiting monodromy around x = oo, the FZZT 
partition function now exhibits what is known as Stokes' phenomenon. As we discussed 
in the introduction, this is the phenomenon whereby the analytic continuation of the 
asymptotics of a function in one region does not correctly reproduce the asymptotics of the 
function in another region. (See also appendix B for a brief review of Stokes' phenomenon, 
summarizing ||39|| .) 

The Airy function is actually a paradigmatic example of Stokes' phenomenon. In the 
region a; -|-oo, the Airy function is given approximately by 
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Attempting to analytically continue the asymptotics counterclockwise around large x to 

X - 

the Airy function on the negative real axis is actually 



oo, one would find Ai{x) ~ e"'"3**^ ^^^^ there. However, the correct asymptotics of 



i.e. it is a linear combination of the two saddles (|4.1|) . 

The reason this happened is that as we varied x from x — +00 to x — —00, we crossed 
a Stokes' line at 

2n 

arg(x) = ±— (4.4) 

Recall from the introduction that Stokes' lines are the places where various saddle-point 
contributions to an integral appear and disappear. In our example, one can see this by 
starting from the negative real axis, where according to ( [4.3| ), both saddles contribute to 
the Airy integral. Upon crossing the Stokes' lines ( [4.4|) , however, the subdominant saddle 
disappears entirely from the asymptotic expansion, until one reaches the positive real axis, 
where the function is given by (|4.2| ). The Stokes' lines occur at precisely the points where 
the disappearing saddle is most subdominant. 

We can also describe the effect of Stokes' lines in a slightly different way: the presence 
of Stokes' lines implies that in some regions, certain saddle points might not contribute 
at all to the integral. Therefore, the naive procedure of just summing over all the saddle 
points is not valid here. A case in point is again the x +00 asymptotics of the Airy 
function (|4.2|) . There we see that the function is dominated by just one saddle. The other 
saddle clearly does not contribute; if it did, it would contribute an exponentially increasing 
contribution to the Airy function at x —>■ +00. 

Let us conclude this section by briefiy summarizing two general lessons we can learn 
from this example. 

1. The classical saddle-point approximation is certainly valid, but in the exact quantum 
answer, we might need to sum over saddles. Because of Stokes' phenomenon, not all 
the saddles necessarily contribute in various asymptotic regions. Even the dominant 
saddle sometimes does not contribute. lil 



The relation between the various sheets of the Riemann surface and the exact answer was 



explored also in [46|, where the proposal seems to be to sum over all saddles. 
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2. The quantum target space (the complex x plane) differs significantly from the classical 
target space (the Riemann surface). The various unphysical sheets of the Riemann 
surface disappear, also because of Stokes' phenomenon. In a wedge around the erst- 
while branch cut, what were classically interpreted as the unphysical sheets lead to 
exponentially small corrections to the exact, quantum answer. 



5. Comments on Other Backgrounds 

In section 2, we extracted the asymptotics of the FZZT partition function ^/'(x, t) at 
large positive x and small h for general {p,q)- Combining this with the fact that t(^{x,t) 
is an entire function of x, we argued that il}{x,t) exhibited Stokes' phenomenon. Thus, 
we expect that our conclusion above about the disappearance of the Riemann surface 
is generally true and that this phenomenon comes about in a way similar to what we 
saw in the (p = 2, g = 1) model. Certain regions in the unphysical sheets lead to small 
nonperturbative corrections to the semiclassical answer in the physical sheet. Meanwhile, 
other regions in the unphysical sheets do not contribute such effects. 

However, a true generalization of our analysis from (2, 1) to other values of (p, q) 
requires clarification of two issues: first, the role of the ZZ branes, which exist for higher 
(p, q) but not for (2, 1), and second, the overall nonperturbative consistency of models with 
higher (p, q) . 

First, let us discuss the ZZ branes. Recall that at the classical level, the number of 
background ZZ branes is measured by the A-periods of the one-form ydx on the Riemann 



surface M.p,q |^ 



ydx = Nih (5.1) 

where Ni is the number of ZZ branes of type i. In the simplest backgrounds without ZZ 
branes the A-periods vanish and the surface degenerates to a genus zero surface. The mod- 



uli of M.p,q fall into two classes ||2^,^. Moduli which preserve the A-periods correspond 
to closed string backgrounds. Moduli which change the values of these periods arise only 
when ZZ branes are added. 

It is clear that this picture must be modified in the exact quantum theory. Even 
without a nonperturbative definition of the theory it is clear that the periods §^ ydx = Nih 
are quantized when ^ 7^ and therefore they cannot change in a continuous fashion by 
varying moduli. But in order to understand the exact nonperturbative theory, we need a 
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definition of tlie tlieory wfiicli goes beyond tfie worldslieet expansion. We take tlie double 
scaled matrix model to be this definition. 

In the matrix model, the ZZ branes represent eigenvalue instantons, corresponding to 
sub-leading saddles of the matrix integral where some of the eigenvalues are away from 
the cut. Thus, unlike the classical theory, which is characterized by fixed values of the 
integers Ni, in the exact theory we must sum over the A^^. This sum over ZZ branes is 
automatically incorporated in the exact, nonperturbative matrix integral. In the end, the 
exact answer is characterized only by the closed string backgrounds. 

The second issue concerns the nonperturbative existence of the [p, q) models. It is 
well known that certain values of (p, q) (e.g. {p, q) = (2, 3), corresponding to pure gravity) 
do not exist nonperturbatively ||47| - |49| , p5| , |36[] . This happens when the double-scaled matrix 
model potential is not bounded from below. We can study this problem in the continuum 
using the FZZT partition function and its relation to the effective potential [p2 ,23,26] 



i;{x,t) ^ e-^^^^^''^/^'' (5.2) 

For large \x\ away from the negative real axis we have from ( |2.46|) 

^Veffix) ^ -r]Cx^ (5.3) 

Therefore, the effective potential is bounded from below only when 

T] = — sign(sin((j7r/p)) < (5-4) 

For example, in the {p = 2, q = 21 — 1) models ( |5.4| ) is satisfied only for I odd, and it is 
never true in the unitary models with {p, q = p + 1). 

Notice that for the nonperturbatively consistent models, i(^{x,t) vanishes as a; — > +oo. 
This reflects the fact that the eigenvalues are not likely to be found there. This also specifies 
boundary conditions for the differential equations ( |2.14| ) satisfied by ?/', and these boundary 
conditions are sufficient to determine uniquely iIj{x, t). Conversely, in the nonperturbatively 
inconsistent models the semiclassical value of il){x,t) diverges as x -|-oo. Thus there is 
a nonperturbative ambiguity in the definition of il){x^t) corresponding to the freedom to 
add small exponential corrections to the dominant contribution ( [5.2| ). So we see how the 
nonperturbative problems of these models, which are associated with the unboundedness 
of the potential, manifest themselves here in the ambiguity of defining the FZZT partition 
function 'i/j{x,t). 
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6. Summary and Discussion 

We have explored the relation between the semiclassical geometry seen by the FZZT 
branes in minimal string theory and the exact results as computed by the matrix model. We 
have seen that the various sheets of the Riemann surface correspond to different saddle 
points of the effective theory on the FZZT brane. For some ranges of the value of the 
boundary cosmological constant all the saddles contribute to the answer, while for some 
other ranges only a subset of all saddles contribute. The precise matrix model definition of 
the theory tells us which saddles contribute and which do not contribute. So the Riemann 
surface, which plays a crucial role in the perturbative analysis of the model, suffers drastic 
modifications when we consider the full nonperturbative aspects of the theory. 

We have also given a quick derivation of the relation between the Kontsevich matrix 
model and the ordinary double scaled matrix model, clarifying the relation between the 
various open string descriptions of the theory. The Kontsevich model arises as the effective 
theory on the FZZT branes after integrating out all the open strings corresponding 
to the ZZ branes. The degrees of freedom of the Konsevich model are, roughly speaking, 
"mesons" made out of the fermionic strings stretched between FZZT and ZZ branes ||2^ . 

We have discussed in detail the simplest (2, 1) model, but we also argued on general 
grounds that the disappearance of the Riemann surface due to Stokes' phenomenon is a 
feature of all of the (p, q) models. 

Our paper was partially motivated by the discussion of physics behind the horizon in 
3% . We have a somewhat similar problem. The Riemann surface is analogous to space- 



time and the second sheet is analogous to the region behind the horizon. The FZZT brane 
observables could be viewed as probes of the spacetime geometry. Note that the parameter 
X = fiB is controlled in the boundary region "outside" the horizon. By analytically con- 
tinuing semiclassical answers we get to explore the second sheet of the Riemann surface, 
much in the same way that the region behind the horizon is explored in |]3^. The black 
hole singularity, where quantities diverge, is somewhat analogous to the x ^ oo region of 
the second sheet, where again the expectation values of the analytically continued FZZT 
branes diverge. In the exact answer, however, this saddle point ceases to contribute before 
its value becomes very large. In both cases, the holographic theory tells us that there 
are no divergences in this region. We expect that further analysis of these simple exactly 
solvable examples might yield interesting general lessons for how to think about quantum 
gravity in higher dimensions. 
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Another motivation for our work was the resemblance, at least at the perturbative 
level, between the minimal string theories and the topological string (see, e.g. the discussion 
of the [p, 1) models in [0). Nonperturbatively, however, the connection is less clear. While 
we lack a generally accepted nonperturbative definition for the topological string, minimal 
string theories have an exact, nonperturbative formulation in terms of the dual matrix 
model. Given the similarities between the two theories, it is natural to suppose that some 
of the lessons from our work might be relevant to the topological string. Let us just briefly 
mention a few. 

First, we have seen how the semiclassical target space (the Riemann surface) is drasti- 
cally modified by nonperturbative effects. In the topological string, the Riemann surface is 
intimately connected with the target space. The Riemann surface is the surface H{x, y) = 
in C^. The Calabi-Yau is given by the equation uv + if(x, y) = in C^. This Calabi-Yau 
is a C-fibration over the complex (x, y) plane, and the discriminant locus of the fibra- 
tion is the Riemann surface H{x, y) = 0. Our results raise the question of whether in a 
proper nonperturbative definition of the topological string the Calabi-Yau might also be 
drastically modified by quantum effects. 

Another striking feature of our analysis is the role of Stokes' phenomenon. Semiclas- 
sical target space is viewed as a saddle-point approximation to some effective theory on 
the brane probe. Nonperturbatively, we must sum over different saddle-points in a pre- 
scribed fashion. As we saw with the Airy function, the result was that some portions of 
the semiclassical target space contributed to physical observables, while others did not. It 
even happened sometimes that the dominant saddle-point did not contribute. It will be 
interesting to see if these phenomena play a role in the nonperturbative topological string. 

A third possible application to the topological string is the role of the ZZ branes. These 
correspond to eigenvalue instantons in the matrix model. The classical vacuum of the 
matrix model corresponds to placing all of the eigenvalues into the dominant minimum of 
the matrix model potential. However, in the exact answer we must sum over all vacua (i.e. 
sum over all instantons), obtained by filling the other critical points of the potential with 
any number of eigenvalues. In the continuum limit, this means that we must integrate over 
a subset of the moduli of the Riemann surface describing the normalizable deformations due 
to ZZ branes. Note that the closed-string couplings are non-normalizable deformations, 
and hence we do not integrate over them. This suggests that in the nonperturbative 
topological string, one should also integrate over some of the moduli of the Calabi-Yau 
(the normalizable modes) but not others. 
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In the context of the topological string, it has been suggested that the Riemann surface 
is covered by patches and the D-branes in different patches are related by (generalized) 
Fourier transform |3^,46|. Comparison with the two matrix model suggests, as mentioned 
in footnote 3, that the theory has two distinct branes det(a; — M) and det{y — M). As 
in [^, these are natural in different patches on M.p,q consisting of the first sheet of the 
x-plane and the first sheet of the y-plane (note that these two patches do not generally 
cover the whole surface). As in [^^, in the classical theory these branes are related by 
Legendre transform, with the boundary cosmological constant x and its dual y satisfying 
the defining equation of M.p,q- It is likely that nonperturbatively they are related by a 
Fourier transform. This can be interpreted as a relation between different branes rather 
than as a relation between different patches of the surface. 

Finally, let us mention a more mathematical potential application of our work. We 
have seen that the asymptotics of the Baker- Akhiezer function ( p.31| ) shows very clearly 
the emergence of the classical Riemann surface through the one-form ydx. It is an old idea 
||36|| that the full Baker- Akhiezer function should be used to define a "quantum Riemann 
surface," associated with the string equations [P, Q] = /i in a way analogous to the associ- 
ation of a Riemann surface to the stationary KdV equations, in which case [P, Q] = 0. A 
closely related point is the relation of the matrix model partition function and KdV flows 
to the inflnite Grassmannian. In particular, in the free fermion interpretation of the infi- 
nite Grassmannian one needs to introduce an operator which does not create monodromy, 
(such as twistfields in conformal field theory) but rather Stokes multipliers. Such opera- 



tors, called "star operators" in [|35| , |3q1 , are not at all well-understood. It was suggested 
in []35|j3^ that the point in the Grassmannian created by star operators should define a 
"theory of free fermions on a quantum Riemann surface." A similar suggestion has recently 



been made in p3]. Perhaps it is a good time to revisit these issues. 
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Appendix A. Geometric Interpretation of the Lax Formalism 

In this appendix, we will study the Lax operators Q and P in the semiclassical h ^ Q 
limit. Much of this section will consist of collecting and streamlining many facts that are 
scattered throughout the literature. However, in the process of organizing this material, 
several new insights will emerge. 

Our main goal is to provide a geometric interpretation for P and Q in terms of the 
Riemann surface M.p^q of minimal string theory. How this geometric interpretation is 
modified at 7i 7^ is an important question. In , it was proposed that by generalizing 
the Burchnall-Chaundy-Krichever theory of KdV flow, phrased in terms of framings of 
line bundles, to framings of a flat holomorphic vector bundle over the space of x,tj, one 
could define a notion of a "quantum Riemann surface." It would be nice to understand 
the relation of this proposal to the geometrical interpretation given below. 

A.l. A brief review of the Lax formalism 

First, let us take a moment to recall briefiy the definition of the Lax operators Q and 
P of minimal string theory. (For a more thorough review, see e.g. []T^.) These operators 
are a convenient way to package neatly the data (physical correlation functions) of minimal 
string theory. They are differential operators, of degree p and q respectively, in 



where r = ti is the coupling to the lowest- dimension operator. Explicitly, we have 



where the coefficients Up-j{t) and Vq-j{t) represent various two-point functions of physical 
closed-string operators. They depend on the closed-string couplings t = (ti,t2, ■ • •)• -^^^ 
instance. 



d — hdr 



(A.l) 



i=2 





Mp_2(t) oc d^logZ 



(A.3) 



corresponds to the "specific heat" of the string theory. 
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To solve minimal (closed) string theory, 
of Q and P on the closed-string couplings t 
Q and P satisfy the string equation 



, we simply need to solve for the dependence 
= (ti, ■ ■ •)• This is done by requiring that 



[P, Q] = h (A.4) 

along with the KdV flows 

h^^^=UQ^/% n|^ = [P,Q^/n (A.5) 
The compatibility of the latter with the former implies that P is given in terms of Q by 

P= + (A.6) 

fc>l 
k^O mod p 



Substituting this back into ( [A .41 ) then gives a set of coupled differential equations for the 



coefficient functions of Q. These equations can be solved order by order in h, resulting in 
a perturbative expansion for P and Q 

Q = Qo{d, t) + hQi{d, t) + Q2{d, t) + ... 

(A.7) 

p = Po{d,t) + hPi{d,t) + P2{d,t) ^ . . . 



where by convention the operators on the RHS of ( A.7 ) are ordered such that the (i's are 
all on the right. 

A. 2. The Lax operators in the semiclassical limit 

Now let us take to obtain a much simpler set of equations for the Lax operators. 
In this limit, the string equation ( [A.4| ) becomes 

dPo dQo dPo dQo 



dd dr dr dd 



1 (A.8) 



i.e. the commutator is replaced with a Poisson bracket. To see this, note that every time 
d = hdr acts on something to its right, it contributes a factor of h. Therefore the leading 



order contribution to the commutator is the Poisson bracket ( |A.8|) . 

The solution to this equation is well-known (see e.g. section 4.5 of [0). It is simply 

Po{d;t) =y{x;t) with x = Qo{d;t) (A.9) 
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where y{x;t) is the singular part of the large matrix model resolvent in the closed-string 
background labelled by t. Since (x, y) lie on the Riemann surface Mp,q, ( |A.9| ) implies that 
aX h = 0, the simultaneous eigenvalues of Qo and Pq also lie on the same Riemann surface, 
i.e. 

iQo:Po)eMp,q (A.IO) 

The fact that they can be written in the form ( [A.2| ) as polynomials in d implies that the 
eigenvalue of d is the uniformizing parameter for Aip^q. Thus we can write ( |A.9| ) as follows: 

Qo = x(2 = rf;t), PQ=y{z = d;t) (A.ll) 

Thus we have reduced the algebraic-differential problem of solving the genus zero string 
equation to the geometric problem of finding the uniformizing parameter of M.p,q- This 



problem has been solved in various special cases. For instance, in [23], it was found that 



x{z)=Tp{z), y{z)=Tq{z) (A.12) 

in the conformal background. (To keep the equations simple in this appendix, we will 
rescale y so as to remove the coefficient C. This will have no effect on arguments below.) 

Although it is in general a nontrivial exercise to extract from the string equation the 
higher order h corrections to the Lax operators, it is actually easy to obtain the first order 
h corrections Qi and Pi. This is because the coefficient functions of Q and P, being 
closed-string observables, have an expansion in ti? (the closed-string coupling), not h (the 
open string coupling). (Note that this statement is only true for the particular ordering 
prescription we used in defining the Lax operators (|A.2|) .) Thus Qi and Pi arise only from 
the non-commutation of d and the coefficient functions. This gives 

Qi = ^drd^x{z,T)\^^^, Pi = ^drd^y{z,T)\^^^ (A.13) 

Here we have used (|A.llj ), and, as noted above, Qi and Pi are defined with the d's all on 
the right. 

Finally, we should note that the discussion of P and Q in this appendix is limited to 
the classical backgrounds without ZZ branes, where the surface A4p^q has genus zero and 
a number of pinched cycles. It will be interesting to see how to generalize this discussion 
to backgrounds with ZZ branes present. Then the pinched cycles of Mp,q are opened 
up and the surface no longer has genus zero. In such backgrounds, z is no longer a 
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good uniformizing parameter, and our interpretation of P and Q will have to be modified 
accordingly. 



KdV flow and deformations of Aip^q 

Having shown that the simultaneous eigenvalues of Q and P (we will drop the subscript 
from this point onwards) are nothing but the coordinates {x,y) of M.p,q-, we can now 
provide a geometric interpretation of the KdV flow equations ([A.5| ). The KdV flows tell us 
how to deform Q and P from a closed-string background t to a nearby background t + 5t. 
This gives rise to a deformation of M.p,q- Therefore, on general grounds, the genus-zero 
KdV flows must be equivalent to the singularity-preserving deformations of M.pn discussed 



m 



We can check our claim explicitly in the conformal background. After a lengthy 
calculation, whose details we will skip, one derives the following deformations of P and Q 
from the KdV flow equations ( |A.5|) : 



dQ 



Up-i{d) 



and 



dP 

dTr,s 



Uq-lid) 



Tqr{d)Ups-l{d) 
Up-l{d)Uq-M 



Tps{d)Uqr-l{d) 
Up-l{d)Uq-l{d) 



Uq-lid) 



Uqr- 



ps- 



lid) 



Uq-l{d)Up-l{d) 



(A.14) 



(A.15) 



in the conformal background, up to an overall normalization factor. Here r^^g is the coupling 
associated to the continuum operator Vr,s; it is related to the matrix model couplings tj 
by a linear transformation. (The change of basis between matrix model and continuum 
couplings is discussed in [^.) It is important that both ( |A.14| ) and ( |A.15[ ) are polynomials 
in d; this is required by the deflnition of Q and P. Note also that the degree of the 
deformation to Q is always less than p, but there is no restriction on the degree of the 
deformation to P. 

Since the curve for M.p,q in the conformal background is 

F(Q, P) = Tq{Q) - Tp{P) = (A.16) 

the deformation to the curve due to (|A.14[) -( [A.15| ) is 



dF 



Uq-l{Q)Up-lid) 



Tps{d)Uqr-l{d) 



+ Up-l{P)Uq-l{d) 



Up-l{d)Uq-l{d) 

Uqr—ps — 1 (^d) 



Up-liP)Uq-l{d) 



Tqr{d)Ups-l{d) 
Up-l{d)Uq-l{d) 



Uq-lid)Up-l{d) 



Upq-l{d) 



Tps{d)Uqr-l{d) - Tqr{d)Ups-l{d) 
Up-l{d)Uq-l{d) 



(A.17) 
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In the second line, we have substituted ( |A.12|) for Q and P and we have used the identity 



Um-i{Tn{z)) = Umn-i{z) /Un-i{z) . Further use of this identity leads to 

U,.i{Q)T,{Q)Ur-i{P) - Up-i{P)Tr{P)Us-i{Q) (A.18) 



which agrees exactly with the singularity-preserving deformations of M.p,q found in |23]. 



This confirms very explicitly the equivalence between the KdV flows and the deformations 

of Mp^q. 

We should mention that for p = 2, the equivalence of the KdV flows and the 
singularity-preserving deformations of M.p,q can be seen more directly using the formulas 
in ||3^,^. There it is shown, using the representation of the KdV equations as flrst-order 
matrix equations, that one can deflne an "^-deformed" Riemann surface = F{x;t,h) 
which reduces as ^ to the classical Riemann surface (what we call A^p,g) of the matrix 
model. Here F{x;t, K) is a polynomial in x, which depends in a complicated way on the 
closed string couplings t. Although we will not discuss the details here, one can show that 
at /i = 0, the Riemann surface reduces to 

y"^ = {x + u{t)){B{x] t)f (A.19) 

where B{x,t) is a polynomial in x as well as in the Gelfand-Dickii potentials Rj[u\- (See 
eq. (2.35) of [^5|.) The form (|A.19|) shows immediately that the KdV flows are singularity- 



preserving deformations of M.p^qi since as we change the couplings t, the RHS of (|A.19| ) 
always has only one branch point at a; = —u{t) and singularities at the roots of B{x;t). 

It is interesting to contrast this with the Burchnall-Chaundy-Krichever theory of sta- 
tionary KdV flows. There the Riemann surface is obtained from simultaneous eigenvalues 
of the differential operators [P, Q] = 0. The KdV flow preserves the Riemann surface 
moduli and instead is straight-line flow along the Jacobian of the Riemann surface ||5^ . 



Instantons and the singularities of Aip^q 

Finally, we will discuss the connection between instantons and the singularities of 
■Mp,q- Instantons were studied using the classical limit of the Lax formalism by Eynard 
and Zinn- Justin in |]5T|]. Let us briefly review the logic of their analysis. To leading order, 
an instanton corresponds to an exponentially-suppressed perturbation e{t) of the speciflc 
heat u{t) and all other physical correlation functions. Thus in the h ^ limit, we can 
write 

e'/e = (A.20) 
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for some constant r which measures the strength of the instanton. (The derivative in 
( |A.20|) is with respect to the lowest-dimension couphng r.) Since as — > we can ignore 
the t dependence of u{t), we might as weU set u{t) = 1. Then ( [A.20|) can be written as 
de = e{d + r), which imphes that 

f{d)e = ef{d + r) (A.21) 

for any function f{d). 

The next step in the analysis of is the observation that the instanton deforms the 



Lax operators by 

5Q^eS{d), 5P = eR{d) (A.22) 

where S{d) and R{d) are polynomials in d of degree p — 2 and q — 2 respectively. Since 
this deformation must preserve the string equation [P, Q] = h, this leads to the following 
constraint at linear order in e: 

[P, SQ] + [5P, Q] = (A.23) 
Substituting ([O^ ) and using (^^), we find 



iP{d + r) - P{d))S{d) = {Q{d + r) - Q{d))R{d) (A.24) 

This constraint must be satisfied for every d and for some constant r. Since Q{d^-r) —Q[d) 
and P{d + r) — P{d) are degree q — 1 and p — 1 respectively, but S{d) and R{d) are only 
degree p — 2 and q — 2 respectively, ( |A.24| ) implies that Q{d-\-r) —Q{d) and P{d + r) —P{d) 
must share a common root. Thus there exists some d = do where 

{Q{do + r),P{do + r)) = {Q{do), P{do)) (A.25) 



The authors of |^ use ( [A.25| ) to solve for r, and then use ( |A.24| ) to solve for S{d) and 
R{d). 

With the geometric interpretation of the previous sections in hand, we can offer some 
new insights into the analysis of the instantons. The condition (|A.25| ) is equivalent to the 
condition that J^p,q have a singularity (pinched cycle) at the point (x, y) = (Q((io), P{do))- 
This shows that the instantons are in one-to-one correspondence with the singularities of 
Aip^q. It confirms in a direct way the analysis of and the interpretation of the ZZ 
branes as instantons. 
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In [g3[ , it was also argued that the period of y dx around the S-cycle passing through 
the (m, n) singularity computes the (m, n) instanton (ZZ brane) action, i.e. 

Zm,n<x. j) ydx= / y{z)x'{z)dz (A.26) 

with the constant of proportionality independent of m and n. The derivative of this with 
respect to the lowest-dimension coupling r must then be essentially the constant r defined 
in ( |A.20|) . Indeed, a calculation similar to (|2.36|) shows that 



dz = rm,n (A.27) 

- m ,n 

as expected. This provides a non-trivial check of the formula (|A.26|) for the instanton 



actions derived in ||2^. It also generalizes (and simplifies) the analysis of [0, where 
( |A.27 ) was proven for the special case of the conformal background. 



In the conformal background, one can check that the instanton actions rm,n are always 
real. However, in a general background they will be complex. For instance, in the (2, 2m — 
1) models perturbed by the lowest-dimension operator, one can use the formulas in 
to prove this explicitly for m odd. When the rm,n are complex, the corresponding (p, q) 
minimal string theory is expected to be nonperturbatively consistent and Borel summable. 
In these cases, the rm,n come in conjugate pairs, so that even though they are complex, 
the total instanton correction to the partition function is real. 



Appendix B. A Brief Review of Stokes' Phenomenon 



In this appendix we will briefly review Stokes' phenomenon, summarizing |3^. Con- 
sider the following integral 

I{x) = [ ds e-*'^^^'") (B.l) 

J Co 

where S{s, x) is holomorphic in s, and Cq is a contour in the complex s-plane, chosen so 
that the integral exists and admits an analytic continuation to some region of the complex 
x-plane. We are interested in the ^ ^ asymptotics. 

Since ^ = 0, lines of constant ImS are perpendicular to lines of constant Re 5; 
i.e. they are gradient lines of KeS. We would like to deform the contour Cq in ( [B . 1| ) 
to a steepest descent contour C - a gradient line of Re 5 along which ImS is constant. 
(The latter requirement prevents cancellation between different non-saddle portions of the 
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contour in the leading ^ — >■ approximation.) At a generic point such hues do not intersect. 
However, the saddle points ^ = are characterized by having two intersecting steepest 
descent lines. 

Since typically the different saddle points occur at different values of ImiS, it is im- 
possible to deform the contour Cq to a steepest descent contour (constant ImiS) C passing 
through all of them. However, if the steepest descent contours Ci^2 through two different 
saddles labelled by 1 and 2 pass near each other, and have the proper asymptotic behavior, 
we can deform Co as follows. We deform it to a steepest descent contour C which starts 
close to Ci passes near the saddle point 1, then passes near the saddle point 2 and finish- 
ing close to €2- Such a contour must be compatible with the asymptotic behavior of the 
original contour Cq. Alternatively, if Ci and C2 asymptote to each other at infinity and 
Re5 ^ +00 there, we can take C = Ci + C2 (see figure 1 and the example below). This 
makes it clear that the two saddles contribute to the integral. 

Now let us vary the parameter x in (|B.1|) and examine the saddles and the contour 
C. There are two interesting things that can happen. The first, more trivial phenomenon 
is when the two saddles exchange dominance. This occurs across lines in the complex 
X plane called "anti-Stokes lines," where the values of Re 5 at the two saddles are the 
same. The second, more interesting critical behavior happens across the "Stokes lines," 
where the values of ImiS at the two saddles are the same and the topology of Ci^2 changes. 
Beyond this point the contour with the correct asymptotic behavior, or equivalently a 
smooth deformation of the previous contour C, does not pass through the two saddles but 
only through one of them. It is possible to find another contour which passes through 
both of them, but it does not have the correct asymptotic behavior. The exchange of 
dominance of two saddles and the abrupt disappearance of the saddle-point contribution 
to the integral I{x) both contribute to Stokes' phenomenon. As mentioned in the body of 
the paper, this is the phenomenon in which the analytic continuation of the asymptotic 
expansion of a function does not agree with the asymptotic expansion of the functions' 
analytic continuation. 
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(A) (B) 

Fig. 1: The steepest descent lines Ci,2 pass through the saddles points 1 and 
2. The dotted line is the original integration contour Cq. For one value of x the 
situation is as in figure A, and Co can be replaced by C = Ci + C2 because the two 
contours C\ and C2 meet at an asymptotic infinity where the integrand vanishes. 
Then the integral receives contributions from the two saddles. For another value 
of X, as in figure B, the steepest descent contour is given by C = C2 alone, and so 
the integral receives a contribution only from the saddle 2. The transition occurs 
for the values of x for which Im5(l) = Im<S(2). 



As an example, consider the Airy function 

/+00 3 
rfse*(^+^«) (B.2) 
-00 

The behavior as |s| 00 allows us to deform the contour to start in the wedge ^ < 
arg(s) < TT and end in the wedge < arg(s) < |-. The two saddles at s ±^/—x are as in 
figure 1. There is an anti-Stokcs line located on the negative x axis. Here the two saddles 
are purely imaginary (i.e. RciS = 0) and they exchange dominance. One can also check 
that the lines | arg(a;)| = ^ are Stokes lines. Thus, Figure A applies to ^ < | arg(a;)| < tt 
and Figure B applies to < | arg(a;)| < As one crosses the Stokes lines starting from 
the negative real axis, the dominant saddle ceases to contribute. 



Appendix C. Numerical Analysis of {p, q) = (2, 5) 

In this appendix we will analyze in detail the example of {p, q) = (2, 5), using numerical 
methods where necessary. The purpose of this analysis is mainly to verify that the lessons 
we learned from the example of (2, 1) indeed carry over to more complicated models. 

To begin, we define the Lax operators to be 

2 

Q = d''-U{T), P = J]t2fc+3Q++'/', d=hdr (C.l) 

fe=0 
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This describes a perturbation around the (2, 5) multi-critical point. Let us set tj = —8/5 
without loss of generality. Then the string equation [P, Q] = h takes the form 

+ ^t^u^ -hu-T- (2w'2 + 4uu" + hu") + Y^^^w^^^ = (C.2) 

The Baker-Akhiezer function is determined by the differential equations 

= ^•0, ^V' = ^^xV' (C.3) 

together with the condition that is real and exponentially decreasing as x ^ +00. 

Before we proceed to solve ( |C.2|) and (|C.3|) numerically, let us first discuss the classical 



limit — > 0. At = 0, P and Q take the form 

Q = (f-u, P = -\(f + {Au + U)d^ -{?,u^ + ^hu-h)d (C.4) 

with u{t) the solution to (|C.2| ) with h = Therefore they lie on the Riemann surface 
described by the algebraic equation 

= ^(Q + [sQ^ -a(u+ Q + ?,u' + h^u - Sts] (C.5) 



Here we see explicitly how the Riemann surface takes the form ( |A.19D for all values of the 



closed-string couplings r^ts.t^. In particular, the Riemann surface always has a branch 
point at Q = —u and singularities at the other roots of the RHS of ( P.5|) . Therefore the 
KdV fiows, which change the values of the closed string couplings, indeed correspond to 
singularity-preserving deformations of the Riemann surface. 

Now let us discuss the numerical solution of the string equation ( |C.2|) and the Baker- 
Akhiezer equations ( |C.3|) . The string equation for perturbations around the (p, q) = (2, 5) 



critical point was solved numerically in P7| , |49|] . Here we will repeat the analysis of [^ , |49|] to 
obtain the specific heat u{t) for various values of h. We will then take the analysis one step 
further by numerically solving ( |C.3|) for the Baker-Akhiezer function. For simplicity, let us 
limit ourselves to the conformal background perturbed by the lowest-dimension operator. 
The conformal background corresponds to 

r = 0, t3 = l, t5 = (C.6) 



Up to a trivial shift of u and r this is identical to the setup considered in ||4^ . To see 
that this is the conformal background, simply substitute ( |C.6| ) into the formula ( |C.5| ) for 
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the Riemann surface. Since the string equation (|C.2|) is solved by u(r = 0) = 1 (modulo a 
discrete choice for the root of the cubic polynomial) the curve becomes 

y^ = ^{x + l)i4x^-2x-lf (C.7) 

which is indeed the same as T2{y) = T^(x) after a rescaling of y. 

Shown in figure 2 is the specific heat u{t) versus r for various values of h. At large |r| 
the specific heat asymptotes to the classical solution Uci{t) ~ sign(r)|r|^/"^. Meanwhile, 
at small |r| the specific heat oscillates faster and faster as h is decreased, since here the 
function is trying increasingly hard to interpolate smoothly between the classical discon- 
tinuity Uc\{t = 0) = ±1 at r = 0. Evidently, the classical limit of u{t) is not well-defined 
for small |r|, although the quantum answer is smooth. 




Fig. 2: The specific heat u{t) as a function of the lowest-dimension couphng r, 
for h = 1,0.5,0.3. These plots were obtained by numerically solving the string 
equation (|C.2|) in the conformal background (C.6). 



Figure 3 contains a plot of the Baker- Akhiezer function i(j{x, t), again for various values 
of h. (The different solutions have been rescaled in order to aid the presentation.) From 
the figure, it is clear that 'i/'(x,t) is decreasing at large positive x, while it is oscillatory 
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for X < ~1. Also, the function is clearly always smooth and real- valued. The bump at 
X ^ in figure 3 corresponds to the location of the (1, 2) ZZ brane, while the trough at 

X ~ is the location of the (1, 1) ZZ brane. As h decreases, the oscillations at a; < — 1 

become faster, and the bump at x ~ becomes more well-defined. This behavior is 

all qualitatively consistent with the leading-order WKB approximation 




-2-1012 



Fig. 3: The Baker- Akhiezer function 'il){x,t) versus the boundary cosmological 
constant x, for h = 1, 0.5, 0.3. The plots have been rescaled for the different values 
of h, so as to improve the presentation. 

A more quantitative comparison between the WKB approximation and the exact 
answer is shown in figure 4, h = 0.3. We see that they are in excellent agreement, except 
for a small region around x = —1 where we expect the WKB approximation to break down 
anyway. 

It should be clear from the discussion that these numerical results confirm many of 
the general arguments in the text regarding the properties of the Baker- Akhiezer function 
ifjixjt). Let us just mention a few. First, "(p^x.i) obviously exhibits Stokes' phenomenon: 
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the analytic continuation of the asymptotics (|C.8|) away from large positive x, where i/j{x, t) 
is exponentially decreasing, leads to the wrong answer for x < —1, where i(;{x,t) is oscil- 
latory. Second, notice that the analytic continuation of the WKB approximation from 
large positive x is accurate up until a; ~ — 1 . The failure of the analytic continuation of the 
WKB approximation beyond x = — 1 is due to the level crossing phenomenon, which results 
in the oscillatory behavior of il){x,t). These facts agree well with the general discussion 
in section 2. Finally, note that the Baker- Akhiezer function is exponentially decreasing 
at large positive x. From section 5, we know that this is the expected behavior for the 
nonperturbatively consistent (2, 5) model. 

3.5 

3 
2.5 

2 
1.5 

1 ' 
0.5 



-0.5 
-1 
-1.5 

-2-10 1 2 

X 

Fig. 4: A comparison of exact Baker- Akhiezer function and its leading-order WKB 
approximation, for h = 0.3. The two are clearly in excellent agreement, except in 
a smah region around x = —1 where the WKB approximation is expected to break 
down. 
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